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Abstract 

Measurements of the magnetoresistivity of graphite with a high degree of control of 
the angle between the sample and magnetic field indicate that the metal-insulator 
transition (MIT), shown to be induced by a magnetic field applied perpendicular 
to the layers, does not appear in parallel field orientation. Furthermore, we show 
that interlayer transport is coherent in less ordered samples and high magnetic 
fields, whereas appears to be incoherent in less disordered samples. Our results 
demonstrate the two-dimensionality of the electron system in ideal graphite samples. 
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Recent experimental and theoretical studies of graphite have renewed the in- 
terest in this system [1,2,3,4,5,6,7,8,9,10]. Experimental results show that, con- 
trary to the common belief, the transport and magnetic properties of graphite 
cannot be accounted for by semiclassical models. Experiment and theory raise 
a number of questions concerning the coupling between the graphite layers. 
The understanding of the transport properties is of primary interest and can 
provide a fundamental contribution to the physics of two-dimensional (2D) 
systems in general. In this letter we deal with two important open questions: 

1. A MIT appears both in the in-plane [1,3,4] and out-of-plane [5] resistivity 
induced by a magnetic field B applied perpendicular to the graphite layers, i.e. 
£?||c-axis. Based on magnetization data [2] and the found 2D scaling similar 
to that in the MIT of Si-MOSFETs (metal-oxide-semiconductor field-effect- 
transitor) [11], Mo-Ge [12] and Bi-films [13], the MIT in graphite has been 
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interpreted in terms of superconducting fluctuations [1,3]. It has been shown 
that the density of states at the Fermi level is enhanced through topological 
disorder, thus leading to the possible occurrence of localized ferromagnetism 
and superconductivity [7]. 

Other interpretation of the MIT in graphite uses the idea that a single graphite 
layer is the physical realization of the relativistic theory of (2+l)-dimensional 
Dirac fermions due to the linear dispersion in the spectrum of quasiparticle 
excitations in the vicinity of the corners of the Brillouin zone [14]. Taking 
this into account and the large Coulomb coupling constant for graphite [5,8], 
a magnetic catalysis (MC) has been proposed [9]. This original explanation 
assumes that a magnetic field perpendicular to the graphite layers breaks the 
chiral symmetry and opens a gap in the spectrum of the quasiparticles at the 
corners of the Brillouin zone. This effect is interpreted as the enhancement of 
the fermion dynamical mass through electron-hole pairing, i.e. a transition to 
an excitonic insulating state. 

Which role does the field direction play? The experimental evidence indi- 
cates that the MIT in Si-MOSFETs occurs independently of the orientation 
of the magnetic field and, therefore, has been considered to be driven solely 
by spin-dependent effects [15]. On the other hand, the MC in graphite would 
be possible only for the case B\\c, i.e. the transition should be absent in the 
parallel case [8]. To our knowledge there is no experimental proof published 
in the literature that the MIT is absent in oriented graphite samples for fields 
applied parallel to the graphite planes. A clear experimental evidence for the 
absence or not of the MIT would give an important hint to search for its origin. 
This is one of the tasks of our experimental work. 

2. In a recent work we have shown that the metalliclike behavior of the out- 
of-plane c— axis resistivity p c is directly correlated to that of the in-plane 
resistivity [5] and that the intrinsic p c {T) of an ideal graphite sample would 
not be metalliclike. These results cast doubts about the interlayer transfer 
integral value of ~ 0.3 eV used all over the literature [20,6] in which neither 
the electron-electron interaction nor charge fluctuations were taken into ac- 
count [21]. In principle, for such a large interlayer transfer one would expect 
coherent transport for the interlayer magnetoresistance at low temperatures. 
In this work we provide an answer to the question whether the c— axis trans- 
port in graphite is or not coherent and whether this (in)coherent transport is 
influenced by sample inhomogeneities. 

We have measured the following samples: a highly oriented pyrolytic graphite 
(HOPG) manufactured at the Research Institute "Graphite" in Moscow and 
denoted as HOPG-3 [3]. The X-ray rocking curve for this sample, see Fig. 
1(a), exhibits a full width at half maximum (FWHM) of (0.64 ± 0.05)°. We 
take this value as a measure for the misalignment of the graphite layers within 
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Fig. 1. (a) X-ray rocking curve of the sample HOPG-3. (b) Angle dependence of the 
in-plane resistivity (•) for the same sample in a magnetic field of 9 T at 2 K. The 
field parallel to the sample at 90° is defined at the resistivity minimum. A second 
run (+) is stopped at 90°. 

the sample with respect to each other. Further two HOPG samples, one from 
Union Carbide Corp. with a FWHM of 0.24° (HOPG-1) and the other from 
Advanced Ceramics Corp. with a FWHM of 0.40° (HOPG-2), have been mea- 
sured. The fourth sample was a Kish graphite sample with a FWHM of 1.6° 
and a in-plane resistivity at 2K and zero field ~ 100 times smaller than for 
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the HOPG samples. AC resistivity measurements were performed by a conven- 
tional four-probe method. The samples were fixed in a rotating sample holder 
inside the bore of a 9T superconducting solenoid. Temperature stability was 
better than 2 mK in the whole temperature range. 

The crucial point of the magnetoresistance measurements in parallel field is 
the misalignment of the sample with respect to the field. The ultimate limit 
for the alignment of the sample would be that of the graphite layers within 
the sample. To adjust the sample to the magnetic field we measured the angle 
dependence of the in-plane resistivity with the current always perpendicular 
to the field. The result is shown in Fig. 1(b) for the HOPG-3 sample. Here 90°, 
i.e. field parallel to the sample, is defined at the minimum of the resistivity. 
The high angle resolution, the strong angle dependence of the resistivity and 
the high sensitivity of its measurement, as well as the excelent reproducibility 
of the absolute angle (of the order of the angle resolution, see Fig. 1(b)) allow 
us to align the sample parallel to the field with an accuracy of ±0.02°, well 
below the FWHM of the rocking curve. 

Figure 2(a) shows the results of the in-plane resistivity with B\\c. We note 
that the metalliclike phase, i.e. dp/dT > 0, is observed between a maximum 
T max and a minimum temperature T m ; n and for fields below a critical field B c 
. The difference between T max and T min decreases as the field approaches B c , 
i.e. T max (i? c ) = T min (B c ), and the metalliclike phase disappears. Figure 2(b) 
shows these data together with fits of T max (B) to the experimentally found [1] 
relation 



where B and B x are free parameters, the first is refered as the offset field [10]. 
The extrapolation of those relations serves to determine the critical field at 
their crossing which for the present case is B^r = (0.12 ± 0.01)T. 

The results obtained for B _L c (aligned as shown in Fig. 1(b)) are shown in 
Fig. 3(a) and the corresponding T min and T max in Fig. 3(b). Whereas T min (B) 
follows a similar relation as in the former case, T max (B) is reproduced by 
Eq. (1) only in the low field range. At high fields the empirical relation 

B 

^max e s 2 is found. We note that the intersection of T max (B) and T min (B) is 
at -Bj = (11.3 ± 0.1)T. The dominant contribution to the misalignment is the 
misalignment of the graphite layers with respect to each other and is taken 
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Fig. 2. (a) In-plane resistivity of the sample HOPG-3 as a function of temperature 
for magnetic fields B = 0, 0.02, 0.04, 0.06, 0.08, 0.09, 0.10, 0.11, 0.12 T (bottom 
to top) applied ||c. Upward (downward ) arrows mark the minima (maxima) of 
the resistivity, (b) Temperatures of the maxima and minima of the resistivity as a 
function of field together with fits for T max : Eq. (1) (-) with B x = 0.13 T; T min : 
dotted line corresponds to Eq. (2) with Bq = 65.9 mT. The intersection of the 
curves gives the critical field B c . 
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Fig. 3. (a) In-plane resistivity of the sample HOPG-3 as a function of temperature 
for magnetic fields B = 0, 1 ... 9 T (bottom to top) applied _L c, aligned as in Fig. 
1(b). (b) Temperatures of the maxima and minima of the in-plane resistivity as 
a function of field together with fits for T m£LX ~. Eq. (1) (-) with Bi = 5.83 T and 
dashed line oc exp(— B/B2) with B2 = 5.11 T; T m - m : dotted line corresponds to 
Eq. (2) (B = 1.59 T). 



to be (0.64 ± 0.05)° (see Fig. 1). The perpendicular component of B\ is then 
(0.13 ± 0.01)T. This value is within the error equal to B^r. From this we con- 
clude that the MIT is, if not solely driven by the magnetic field perpendicular 
to the graphite layers, by far dominated by it. 
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We would like to note two details of the field-driven transition shown in 
Fig. 2(a) for sample HOPG-3, as well as in previous publications [3,4,5]. First, 
in the "insulating" side of the transition at B > B c we would expect that the 
resistivity p — > oo for T — > 0. Instead we observe always a saturation or, upon 
applied field, a weak logarithmic increase of the resistivity decreasing temper- 
ature [16]. The saturation for T — > in the insulating side of the transition 
has been reported for various 2D systems as, for example, in Mo-Ge films [12], 
GaAs/AlGaAs heterostructures [17] and Josephson junction arrays [18]. There 
is no consent on the origin of this saturation. Whatever the reason for it is, 
the main result of this part of the work, i.e. that only the normal component 
of the applied field to the planes drives the transition, remains untouched. 

Second, one would tend to underestimate the real magnitude of the effect 
driven by the applied field because of the relative small change of the re- 
sistivity shown in Fig. 2(a). However, we stress that the relative change in 
the resistivity with field depends on the sample characteristics. Experimen- 
tal data from different graphite samples show a change at low temperatures 
between ~ 20% up to more than one order of magnitude for a field of the 
order of 1 kOe [3,4,5,19]. Qualitatively speaking, the transition is similar for 
all samples studied. 

We discuss now shortly the origin of the MIT in graphite. First we note that 
Eq. (72) from Ref. [10] (T C (B) oc (1 - (BI/B 2 ))VB) fits the data for T min (B) 
as well as Eq. (2) with similar B . In Ref. [10] the authors argue that the 
offset field B ~ nnc/\e\ (n is the charge density) is model independent and is 
related to the minimum field required to fill the lowest Landau level, necessary 
condition to deblock the electrons for pairing and to produce the excitonic gap. 
For the majority carriers n ~ 10 11 cm -2 and B ~ 2 T in clear disagreement 
with the experimental result ~ 0.06 T. But, if we take the minority carrier 
density n ~ 10 9 cm~ 2 [20] we get roughly the measured B . Nevertheless and 
since MC is a 2D phenomenon, it is unclear whether this assumption is valid. 
We note that recent results [5] including those from this work indicate that for 
HOPG samples, the coupling between planes is much weaker than previously 
assumed, casting doubts about the correctness of a 3D Fermi surface with two 
types of carriers with different densities and effective masses for ideal graphite 
[20]. 

Returning to the superconducting scenario which is supported by the magneti- 
zation data [2], we may relate T max (B) to the critical temperature of a system 
of superconducting islands in a semiconducting matrix [1]. In this case meso- 
scopic effects play a role. Indeed, the observed behavior in the parallel case 
shows a change of curvature of T max (B) similar to that found theoretically for 
disordered 2D superconductors [22]. We note also that an exponential decay 
with field for T C (B) has been predicted for 2D superconductors with weakly 
Josephson-coupled local superconducting islands [23]. Because in the parallel 
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Fig. 4. Angle dependence around 90° of the out-of-plane resistivity of the Kish 
graphite (a), HOPG-2 (b) and HOPG-1 (c) samples at B = 9 T and at 2 K. 

case the MIT occurs at large fields, we are able to observe the anomalous decay 
of T max in contrast to the transverse case, see Fig. 2(b). On the other hand we 
may argue that the reason for the difference in the behavior of T max (5) may 
be due to its sensitivity to the intrinsic misalignment of the graphene planes 
of the sample that shows a gauss-like distribution. 

We discuss now the interlayer transport. One possible way to test coherent 
transport across the graphite layers is given by the measurement of a maxi- 
mum in the angle dependence of p c {9) at magnetic fields parallel to the layers 
[24,25,26]. This peak should be absent for incoherent interlayer transport but 
observed if the inequality uj c t > 1 holds, where u c is the cyclotron frequency 
and r the relaxation time of the carriers. Coherent transport means there- 
fore that band states extend over many layers and a 3D Fermi surface can be 
defined. In the other case, incoherent transport is diffusive and neither a 3D 
Fermi surface nor the Bloch-Boltzmann transport theory is applicable [24] . In 
order to check this and the role played by disorder we have performed mea- 
surements of the out-of-plane electrical resistivities with high angle resolution. 



Figure 4 shows the results for three samples at 9 T and 2 K. We observe that a 
weak coherent peak in p c around the parallel orientation (90°) occurs and this 
is larger the larger the FWHM of the corresponding rocking curve. The asym- 
metry seen in the angle dependence is in part due to the small experimental 
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misalignment of the surface of the sample and to the lack of crystal perfec- 
tion. The coherent peak decreases as expected with field. Our results indicate 
that lattice defects not only affect the transport as scattering centers but they 
contribute to enhance the coupling between the layers giving rise to a 3D-like 
electronic spectrum and coherent transport. The absence of coherent peak in 
ideal samples may be related either to incoherent transport or that uo c t < 1 
holds. Although the validity of semiclassical criteria for incoherent transport 
is under discussion [25] we use the Ioffe-Regel-Mott maximum metallic resis- 
tivity p max criterion to evaluate coherent transport [27]. We obtain that only 
for the HOPG-1 and -2 samples p c > p max in the whole T— and B— range, in 
agreement with the absence of coherent peak. 

From our results we draw the following conclusions: The MIT in HOPG is trig- 
gered only by a magnetic field perpendicular to the graphite layers. Therefore, 
it is unlikely that spin effects play a significant role in the MIT. The absence 
of the MIT in the parallel field orientation supports the theoretical approach 
of Refs. [8,9,10], but there is apparently no quantitative agreement [10]. On 
the other hand, the influence of possible superconducting fluctuations on the 
MIT cannot be ignored. The transport perpendicular to the graphite layers in 
highly oriented and less disordered samples appears to be incoherent, demon- 
strating the quasi-2D character of the electron system of graphite. Samples 
defects lead to a better coupling between the layers, a 3D-like behavior and 
coherent interlayer transport, added to the possible local enhancement of the 
density of states which may generate local superconductivity and ferromag- 
netism [7]. 
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